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Abstract 

m 

Consider the nonlinear matrix equation X- £ A*X p, Aj = Q with p, > 0. Sufficient and necessary 

i=l 

conditions for the existence of positive definite solutions to the equation with p, > are derived. 
Two perturbation bounds for the unique solution to the equation with < p, < 1 are evaluated. 
The backward error of an approximate solution for the unique solution to the equation with 
< pi < 1 is given. Explicit expressions of the condition number for the equation with < p, < 1 
are obtained. The theoretical results are illustrated by numerical examples. 

Keywords: nonlinear matrix equation, positive definite solution, perturbation bound, backward 
error, condition number 



1. Introduction 

In this paper the nonlinear matrix equation 



X~Y j A*X p >A i = Q 



(1.1) 



1=1 



is investigated, where A\,A2, .. .,A m are n x n complex matrices, m is a positive integer, pi > 
(z = 1, 2, • • -m) and Q is a positive definite matrix. Here, A* denotes the conjugate transpose of 
the matrix A,- . 

When m > 1, Eq. dl.U is recognized as playing an important role in solving a system of 
linear equations. For example, in many physical calculations, one must solve the system of 
linear equation 

Mx = f, 

where 



/ 





M : 







a: 



a; 






/ 
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arises in a finite difference approximation to an elliptic partial differential equation (for more 
information, refer to (4j )■ We can rewrite M as M = M + D, where 



( x-p< 

o x~p* 

M 



A* A* 

M can be factored as 
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-/ 
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if and only if X is a solution of equation X - 2 A*X Pi Ai — Q. When m — 1 , this type of nonlinear 

i=\ 

matrix equations arises in ladder networks, dynamic programming, control theory, stochastic 
filtering, statistics and so forth 

For the similar equations X ± A'X^A = Q, X s + A*X-'A = QmdX+f l A;X" 1 A,- = / , there 

were many contributions in the literature to the theory, applications and numerical solutions J5[0, 
H [HM IHHil HI [HHH [H [37^3911 . Jia and Gao H derived two perturbation estimates 
for the solution of the equation X-A*X q A = Q with < q < 1 . In addition, Duan et al. J7|] proved 

m 

that the equation X - 2 A*X Si Ai — Q (0 < |(5,| < 1) has a unique positive definite solution. They 

also proposed an iterative method for obtaining the unique positive definite solution. However, 
to our best knowledge, there has been no perturbation analysis for Eq. dl . 11 > with m > 1 in the 
known literatures. 

The rest of the paper is organized as follows. In Section 2, some preliminary lemmas are 
given. In Section 3 , sufficient and necessary conditions for Eq. ( 11.11 ) existing positive definite 
solutions are derived. In Section 4 , two perturbation bounds for the unique solution to Eq. dl. U 
with < pi < 1 are given. Furthermore, in Section 5, we obtain the backward error of an 
approximate solution for Eq. dl.lt with < pi < 1. In Section 6, we also discuss the condition 
number of the unique solution to Eq. dl.lt . Finally, several numerical examples are presented in 
Section 7. 

We denote by C nx " the set of n x n complex matrices, by < J-(" X " the set of n x n Hermitian 
matrices, by / the identity matrix, by i the imaginary unit, by || ■ || the spectral norm, by || • \\p 
the Frobenius norm and by A m - dx (M) and A m i n (M) the maximal and minimal eigenvalues of M, 
respectively. For A = (a\, . . . , a n ) = (ajj) £ C" x " and a matrix B, A ® B — (a/yB) is a Kronecker 
product, and vecA is a vector defined by vecA = (a[, . . . , a^) T . For X,Y e < H nx " 7 we write X > Y 
(resp. X > Y) if X - Y is Hermitian positive semi-definite (resp. definite). 



2. Preliminaries 



Lemma 2.1. Md/- If A > B > and < y < 1, then A^ > B y . 
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Lemma 2.2. [20 ]. For any Hermitian positive definite matrix X and Hermitian matrix AX, we 
have 



(i) X q = — X*(AI + Xy x X'^A q ~ l dk, < <? < 1; 

n Jo 



/-*oo 

(a) xi = smqn I xhai + xy l x(Ai + xy l xh q - l dA, 0<q<l. 

(1 -q)n Jo 
In addition, ifX + AX > (l/v)X > 0, then 
(Hi) \\X-^A\{X + AX) q -X c ')AX-'i\\ < (1 -q)(\\X-iAXX-i\\ + vWX-iAXX-ifW^-AX-iW 2 . 

I — - m 

Lemma 2.3. [7]. The matrix equation X — £ A*X s 'Ai — Q (0 < |<5,| < 1) always has a unique 

i=l 

positive definite solution X. The matrix sequence X^ : 

m 

X s+m +i = Q + 2 A i^A' * = 0,1,2,-.., (2.1) 

converges to the unique positive definite solution X for arbitrary initial positive definite matrices 
X\,X2, ■ ■ ■ ,X m . 

3. Positive definite solutions of the matrix Eq.(l.l) 

In this section, sufficient and necessary conditions for the existence of positive definite solu- 
tions of Eg. ( II. Il l are obtained. 



Theorem 3.1. Eq. ALU has a positive definite solution X if and only if the coefficient matrices A, 
can be factored as 

Aj = (W*Wy Pl/2 Yj(W*W) 1/2 , (3.1) 



where W, Yj (i — 1,2, •■ • ,m) are nonsingular matrices, Z — Qi{W*W) z and 



orthonormal. In this case X — W*W is a solution of Eg. ( 17771 ). 



Z ^ 

Y m j 



is column 



Proof. If Eq. dTTTT i has a positive definite solution X, then there exists a nonsingular matrix W, s.t. 
X = W*W. Furthermore Eq. dl.ll ) can be rewritten as 



which implies that 



WW -Y u A*(W*W) p 'A i = Q, 



(w*wy^Q(w*wy^ + ^(W*w)-^a*(W*w)^(w*w)^ A^wwy^ = i. (3.2) 
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Let Y t = (WW) i Ai(W*W)-i, then A t = (W*W)^ L Y i (WW)2 . Moreover, by Z = Q--(WW)—2, 
(13.21 ) turns into 



Z 



z ^ 

Yi 

Y m i 



= L 



which means that 



Z ^ 
Y\ 

Y m ) 



is column orthonormal. 



Conversely, suppose that A,- has the decomposition ( 13.lt . Let X = W*W. Then 

m m 

x-YA-XP'At = w*w- y(w*W)5y;(w*wr^(w*W) Pi (w*w)"^Fi(w*w)5 

!=1 1=1 

m 

= W*W- YiW'Wjh'YiiW'W)! 

i=l 

m 

= (w*w)Hi - y f;f;-)(w*w)5 = (w*w)?z*z(w*wo' 
= (w*w)Hw*wyiQiQHw*wyHw*w^- =q, 

that is X — W*W being a positive definite solution of Eq. dTTTt . 



□ 



Theorem 3.2. If A\, Ai, ■ ■ ■ ,A m are invertible and Q € H, then Eq. ( Ii.il ) has a positive definite 
solution X if and only if Ai can be factored as 



At = {WMUY'ViNU, 



where N > 0, U is unitary and M > UQU* is diagonal, M — N = UQU* and 



orthonormal. In this case, X—U*MU is a solution of Ea. fiL]\ . 



v 2 
v m ) 



(3.3) 



is column 



Proof. If Eq. dTTTI ) has a positive definite solution X, then X can be factored as X = U*MU, where 
U is unitary and M is diagonal. Therefore Eq. dl.lt can be rewritten as 



U*MU - A*(U*MU) p Ai = Q, 



which implies that 



M - UQU* = j[] UA*(U*MU) p AiU* 
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and 



Y(M - UQU*y?UA*{U*MU) p 'A i U t (M - UQU'ys = I. 



(3.4) 



Let N = (M- UQU*)i and V,- = {U* MU)~ A t U*(M - UQU*Y^. Then A, = (WMUy'ViNU 

Vi ^ 



and M - N 2 = UQU*. Eq.tO turns into j V*V; = /, which means that 



y 2 



is column 



orthonormal. 

Conversely, suppose that A, has the decomposition (13.31 1. Let X = U*MU. Then 

m m 

X-J] A *i XP ' A i = U'MU - Y U*N*V*(U*MUy^(U*MU) Pl (U*MUT^ V/NU 

i=l i=l 

m 

= U'MU -Y^U*N*V*ViNU 

1=1 

= U*(M - N 2 )U = Q, 
that is X = U*MU is a solution of Eg. ( 11. 11 1. 
Theorem 3.3. IfX is a solution ofEq. UJ) with < p t < 1, then 



□ 



X > 



^min(G) + ^ ^min(A*A,) /i mi n (2) 



7 =£7. 



Proof. By Lemma |231 Eq. dl.lt with < pi < 1 always has a unique positive definite solution 
X. Then X > 0, it follows that X« > 0. Therefore X > Q. By Lemma|27T]and Eq.(0», we have 

□ 



X > G + Z A*e"A,- > Un in (0 + 2 A^A^iQ) 7 =j87 



4. Perturbation bounds of Eq. dTTB with < p, < 1 

Here the perturbed equation 

m 

Z-^A^A^G, 0</>;<l, (4.1) 
i=i 

is considered, where A, and Q are small perturbations of A, and Q in Eq. dl.lt . respectively. We 
assume that X and X are the solutions of Eq. dl.lt and Eg. (14. It . respectively. Let AX = X - X, 
AG = Q - Q and AA, = A, - A,-. 

In this section two perturbation bounds for the solution of Eq. dl.lt with < pi < 1 are 
developed. The relative perturbation bound in Theorem 14. 1 1 does not depend any knowledge of 
the actual solution X of Eq. dl.lt . Furthermore, a sharper perturbation bound in Theorem l4.2l is 
derived. 

The next theorem generalizes Theorem 4 in i2(ill with m = 1, ||AGII = to arbitrary integer 
m > 1, ||AGI| > 0. 
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Theorem 4.1. Let b = p + \\AQ\\ - 2(1 - Pi W>\\Ai\\\ *= E/? P 'IIAA ; || (2||A,|| + ||AA ; ||). // 

i=l i=l 

< b < 2(J3 - s) and b 1 - 4(/3 - s){s + \\AQ\\) > 0, (4.2) 

then 

< £> £ l|AA,|| + MIAQII = fi, (4-3) 

w/zere 

2s 2 



2 ||AAj||(ft + V^ 2 -MP- s)(s + jjAejj)) * + ^b 2 -4(J3-s)(s + \\AQ\\) 

Proof. Let 

m 

Q = {AX € -H" x " : ||X" 1/2 AXX" 1/2 || < e ^ ||AA,-|| + w||Ag|| }• 

i=i 

Obviously, Q is a nonempty bounded convex closed set. Let 

f(AX) = Y^(A*(X + AX) Pi Aj - A*X Pi Ai) + AQ, AX e £1 

i=l 

Evidently, / :!)h <j-("x>< is continuous. We will prove that /(£2) c £2. 

m 

For every AX e Q, it follows ||X- 1/2 AO~ 1/2 || < g £ l|AA,|| + w||AQ||. Thus 

m m 

(g £ ||AA,|| + a>\\AQ\\)I > X- 1/2 AXX- 1/2 > (-g £ ||AA,|| - <o\\AQ\\)I, 

1=1 i=l 
m m 

(l+g^j MAAiH + wAfiHIDX >X + AX>(l-£,^ HAAill - wAglUDX 

i=l i=l 

According to ( I4.2l i and ( I4.3l l, we have 

e Z i*u + <*u - , ; 2l||Aal + " s «^ ^ < i. 

ti ^ + 2 -4(/?-s)(s + ||Aei|) b 2(fi-s) 

Therefore 

m 

(l- e ^||AA,||- W ||Aeil)X>0. 

1=1 



From Lemma [Z2l and Theorem l3.3l it follows that 

m 

YjA* ((X + AX) Pi - X*)Ai 

WX-iAXX^-W 2 



X~2 



1=1 



\\X---AXX—2\\ + 



l-e EIIAAfll-MIAfill 

1=1 



(^(l-pdWX^AiX-h 2 ) (4.4) 



\\X-iAXX-i\ 



\\X-?AXX-?\\ 2 

m 

1-6 EIIAAiH-MIAfill 



Therefore 



X"/(AX)X"|| 

m 

X~* Y j A*(X+ AX) p <A~i -A*X Pl Ai Xi +XiAQX~^ 

- i=l 

m 

V X-±A*((X + AX)'" - X p )AiX-± + AgjH \\ 
l 

m 

V X" ? [AA*(X + AX) Pi (A; + AAj) + A*(X + AX) Pi AA,] X^ 



IX^AXX^H + 



\\X-iAXX--- 



h\\2 



1-eIIIAAiii-wiiAen 

i=i 



I -Pi 



\\Ai\Y 



■S 

i=i 



||AA,-||(2||A,-|| + HAA.-II 



■(i+e]T||AA ; || + w ||Aeil) 



IIAfill 



i=l 



* (^T-^)(|^H + i <i+fo+ 



ljAg| 



That is f{Q) Q By Brouwer's fixed point theorem, there exists a AX e Q such that /(AX) = 
AX. Moreover, by Lemma 1231 we know that X and X are the unique solutions to Eq. dl.U and 
Eg. (14.11 1. respectively. Then 



||X-X|| ||AX|| \\X l/2 (X- l/2 AXX- l/2 )X 



l/2|| 



\\X\\ 



\\X\\ \\X\\ 

m 

< ||X- I/2 AXX- 1/2 || < q J] ||AA,|| + a>\\AQ\\. 



□ 



Remark 4.1. With 

m 

2(Z HMjIKM + HAAdD + IIAQII) 

q V IIAA/H + wllAgH = — 

6 + V^ 2 -4CS-^ + ||Aei|) 

we gef p 2 ||AA,|| + co\\AQ\\ -> 0/or AQ -> 0, ||AA,|| -» (i = 1, 2, • • • , m). Therefore Eq.fiTJ} is 

;=1 

well-posed. 

Next, a sharper perturbation estimate is derived. 
Subtracting ( 11.11 ) from ( 14.11 ) we have 

m • pea 

AX + y Smp/?r I [(i/ + X)" 1 xU,]*AX[(i/ + Xr I xU / ]^'" 1 ,i/l = £ + /!(AX), (4.5) 
* Jo 

where 

B, = X''A ; , / = 1,2, - ■• ,m, 

m m 

£ = Y (B*AA; + AA-B,) + y AA-Z^AAj + AQ, 

!=1 1=1 

m 

fc(AX) = y [ATAXIA* - A;V,(AX)Mj - AA* V,(AX)A,] , (4.6) 

!=1 

Zi(AX) = smpi7T I X3(^7 + xy l AX(AI + X + AXy l AX(AI + Xy l X^A Pi ~ l dA, 
n Jo 

/-»CO 

7;(AX) = — — | X5(,l/ + xy l AX(AI + X + AXy l xh Pi ~ l dA. 
n Jo 

m IU.II 2 

Lemma 4.1. Let £ < 1. 7%e« f/ze linear operator I, : H nxn -> < H nX " defined by 

Til ■ y-»CO 

LW = W+ y Sm/7f7r I [(/l/ + X)- 1 X'iA i ]*W[(/l/ + X)- 1 XiAi]A p ~ 1 dA > W e •K" x ". (4.7) 
* Jo 

is invertible. 

Proof. It suffices to show that the following equation 

has a unique solution for every V e <J-{ nxn , Define the operator M : < H nxn — > <j-{ nxn by 



MZ = y — ^ J X" l ^iXH^ + X)" 1 X i ZZ5(/l/ + X)- i ZiA / X-5A Pi - 1 cLl ) ZeW" 

Let F = X" 1/2 WX- 1/2 . Thus © is equivalent to 

Y + MY = X~ 1/2 VX~ I/2 . 
8 



According to Lemma l22l we have 



\\MY\\ < V I X-iA*XHAI + Xy l X(AI + X)- l X±AiX-h P! - l dA 

id 71 Jo 



Jo 

l|2 



imi 



< Yjd- Pi )\\x*Atx-* \\ 2 \\y\\ < jr ^iini < urn, 



which implies that ||M|| < 1 and / + M is invertible. Therefore, the operator L is invertible. □ 
Furthermore, we define operators P,- : C" x " — > "K" x " by 

P,Z, = L-\B*Zj + Z*Bi), Zi e C" x ", i = 1,2, ■ • • ,m. 
Thus,we can rewrite (14.5b as 

m m 

AX = L~'Ag + P/AA,- + L _1 (^ AA*X'"AA,) + L _1 (/j(AZ)). (4.8) 

Define 

[|L _1 H= max HL-'WH, ||Pi|| = max ||PiZ||, i = 1,2, - • , m. 
W 6 < H" X " Z e C" x " 

WW = l [|Z|| = l 

Now we denote 

I = HL-'ir 1 , ^ = \\X~\ ti = \\XP% n i = \\P i \\, e=^^W 2 , * = l,2,---,m, 

' i=\ 

e = y||Aei| + 2(«,||AA,|| + |||AA,|| 2 ), a = y ^ 6(2||A,-|| + ||AA/||)||AA,-||. (4.9) 

i= 1 i= 1 

Theorem 4.2. If 

(1 -o-) 2 

cr < 1 and e< — . (4.10) 

£ + cr£ + 29 + 2 + 6»)(cr£ + 0)) 

f/;en 

2e 

ll*-X||< 



1 + e£ - cr + V(l +^e-cr) 2 -4e(£ + 6) 
Proof. Let 

m m 

/(AX) = L" 1 Ag + V P,AAj + L _1 (V AA*X Pi AAj) + L _1 (ft(AJ0). (4.11) 

(=1 !=1 

Obviously, / : f{ nxn —> <j-(" x " j s continuous. The condition ( 14.10b ensures that the quadratic 
equation (£ + 8)x 2 - (1 + - cr)x + e — with respect to the variable x has two positive real 
roots. The smaller one is 

2e 

v - 



1 + e£ - cr + V(l + - cr) 2 - 4e(£ + (?) 
9 * 



Define Q = {AX e <H nx " : \\AX\\ < v}. Then for any AX e Q, by (l4~T0l we have 

2e 



1 + - cr 



= ; fr + o--l | -2(l-<7-)(ftr + g) 

1 + - cr ~ (£ + of + 26»)(1 + £ e - <f) 

It follows that / - X~ l AX is nonsingular and 

1 1 



-X~ l AX\\ < 



l-\\X-^AX\\ ~ l-fl|AX|| 
Using ( 14.6b and Lemma [2T2l we have 

\\Zi{AX)\\ < {\- Pi )\\AX\hX-Y\\{I + X~ l m~ l \\\\X p '\\<^ 2 " AA "" 



||V,(AX)|| < l^'IIIIAXIIIIX-'HIK/ + x-'Axr'u < 



1 - 4TIIAXH' 
[|AX|| 



l-f||AX|r 

m 

P(AX)|| < J] (llA/lhlZ/CAJQH + (2||A,|| + ||AA,||)||AA,||||V,(AX)||) 

i=i 

\\2 



vW , , ||AX|| 2 |IAX|| \ 

Noting d4?9l > and ( 14.111 ). it follows that 

||/(AX)|| < y||Afi|| + YjjuWAAiW + ^||AA,|| 2 ) + y||/z(AX)|| 



gjlAgj 0||AX|| 2 
< e H h 



WI|AX|| WHAXH 

(TV 0V 2 

< e + H = v, 

1 - f V 1 - f V 

for AX e £2. That is /(£2) Q Q. According to Schauder fixed point theorem, there exists AX, e D, 
such that /(AX.) = X*. It follows that X + AX, is a Hermitian solution of Eq J4.ll ). By Lemma 
231 we know that the solution of Eq.gTQ) is unique. Then AX, = X - X and ||X - X|| < £ 3 . □ 



Remark 4.2. from Theorem \4.2\ we get the first order perturbation bound for the solution as 
follows: 

, m 

\\X - M < yllAgll + ^ «,||AA,|| + O (\\(AAx, AA 2 , • • ■ , AA m , A0|| 2 ) , 

(=1 

05 (AAi, AA 2 , ■ • • , AA m , AQ) -» 0. 



Combining this with A4.8i gives 

m 

AX = L l AQ + L 1 ^(B-AAj + AA*B { ) + (9 (||(AAi, AA 2 , ■ • • , AA„„ A0|| 2 ) . 
as (AAi, AA 2 , • • • , AA m , A<2) — > 0. 
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5. Backward error of Eq. dl.l) with < /?, < 1 

In this section, a backward error of an approximate solution for the unique solution to Eq. 
(11. U with < pi < 1 is developed . 

m 

Theorem 5.1. Let X > be an approximation to the solution X of Eq.fiLl). IfT, = 2(1 — 

(=i 

—It — ™ — — 

Pi)\\X~ A,-X~ S II 2 < 1 and f«e residual R(X) = Q + £ A*X P 'A,- - X satisfies 

i=i 

\\R(X)\\ < — fl-mind,^}, "tore 0, = 1 + p- 1 |||TO||-E> 0, (5.1) 

f/ze« _ 

2||X||p- I || 



||X-X|| < //||fl(X)||, w/iere p = 



0i + A /^-4||X- 1 |||M|| 
Proof. Let 

*F = {AX e <H nxn :|| X- I/2 AXX~ 1/2 ||< 8 2 \\R(X)\\}, 
where Q 2 = Obviously, ¥ is a nonempty bounded convex closed set. Let 

I kM I 

m 

g(AX) = J] A* [(X + AX)"' - X»] A,- + #(X). 
i=l 

Evidently g : *P h-> •7Y" x " is continuous. We will prove that gQ¥) c *F. For every AX e T, we 
have 

||X" 1/2 AXX" 1/2 || < 2 ||^(X)||. 

Hence 

X- I/2 AXX~ 1/2 > -6 2 \\R(X)\\I, 

that is 

X + AX > (1 - e 2 \\R(x)\\)x. 

Using ( 15. U . one sees that 

fe iTOn = ^£B^ < 2p^ip(5ii < 1 



0i + ^-^ix-'illTOH 



Therefore, (1 - 6» 2 ||/?(X)||)X > 0. 
According to ( 14.4b . we obtain 



X-^(AX)X-3 || 

||X-iAXX-5|H 



||X-iAXX _ i|| + 



1-0 2 \\R(X)\\ 



i-e 2 \\R(X)\\ 

= «X)||. 
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2 + \\XSR(X)X- 



z + \\x- l \wVQ\\ 



By Brouwer fixed point theorem, there exists a AX e ¥ such that g(AX) = AX. Hence X + AX 
is a solution of Eq. dl.lt . Moreover, by Lemma l2~3l we know that the solution X of Eq. dl.lt is 
unique. Then 

= HAxii < \\x\\\\x-iAxx-k = e 2 \\X\\\\R(X)\\. 

□ 

6. Condition number 



In this section, we apply the theory of condition number developed by Rice H27I1 to study 
condition numbers of the unique solution to Eq. Jl.lt with < pi < 1. 

6.1. The complex case 

Suppose that X and X are the solutions of the matrix equations dl.lt and d4.lt . respectively. 
Let AA = A - A, AQ = Q - Q and AX = X - X. Using Theoreml4~2land RemarkgJ] we have 



AX = X - X = L- 1 AQ + L- 1 ]T(B- AA; + AA-B,) + O (\\(AA U AA 2 , ■ • • , AA„„ AQ)\\ 2 F ) , (6.1) 

!=1 

as (AAi, AA 2 , • • • ,AA m ,Ag) -> 0. 

By the theory of condition number developed by Rice 112711 . we define the condition number 
of the Hermitian positive definite solution X to dl.lt by 



c(X) = hm sup , (6.2) 

1 71 '72 7m P - 



where p and r)u i = 1, 2, ■ • • , m are positive parameters. Taking £ = 77, = p = 1 in (16.2b gives 
the absolute condition number c a t, s (X), and taking £ = , 77,- = \\A\\f and p = ||2IIf m ( 16.2b 
gives the relative condition number c re i(X). 
Substituting d6.lt into ( 16.2b . we get 

m 

IIL- 1 (Ag+ 2(b;aa, + aa;b,))|| f 

c(X) = - max — — — — — 

€ (AAi , A^, . . . , M*, Afi } ) Mi, .... M)|| F 

AAi £ C" x ", Ag € 7Y" X " 

||L- 1 (pff+2i7iW- + 
- max — 

^ (£i,£ 2 ,-- - ,E m ,H)?Q 11(^1,^2, ■•• ,£m,-ff)lb 

EieC nXn ,H e>H nXn 

Let L be the matrix representation of the linear operator L. Then it is easy to see that 

ZSin P/7T I 1 1 r 1 1 . „ 1 

— — [(AI +X)- 1 X5A i ] T ® [(AI + X)- 1 XiA i ]*A Pl - 1 dA. 

,=1 71 J o 
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Let 



ZT 1 = S + i£, 

L _1 <7 ® B*) = L _1 (7 ® (Z Pi Aj)*) = C/,1 + i£2 tt , 
ZT^fif ® 7)n = L~ l {{X p 'Ai) T ® 7)n = t/ /2 + iO !7 , 



5 -E 
2 5 



vecH = x + iy, vecE, = a,- + ibi. 



On + O i2 t/ii - £/ a 



i = 1,2, 



,aLb T m f, M = (E U E 2 , 



(6.3) 



where x,y, a,-, € S, S, t/,i, C/q, fin, £2,-2 e ^"~ x "~ ; / — 1,2, •• • , m, II is the vec-permutation 
matrix, such that 



vec A — II vec A. 



Then we obtain that 



1 WL-'ipH+ZliiB^Ei + E'BMF 

c(X) = - max — 

£ M^O \\(Ei,E 2 ,- ■ ■ ,E m ,H)\\ F 

m 

WpL^vecH + 2 rjiL^dl ® B*.)vecEi + (Bj ® 7)vec£?)|| 

1 !=1 

! M ±0 IKvecEi, vec£ 2 , • • • , vec£ m , vecff)|| 

||p(S + + iy) + 2 i/iKI/a + IQaXoj + + (tfa + IQo)(fl« - tt/) 

- max 

£, m + ll(vec£i , vec£ 2 , • • • , vecE,,,, vecH)|| 

1 \\(p ScijiUumUi,- ■ ■ ,rj m U m )g\\ 

- max 

tg*o WsW 
^ II (pS c , 771C/1, mUz,--- >rimU„ 



EieC nXn ,H e-H" 



Then we have the following theorem. 
Theorem 6.1. 77ie condition number c(X) defined by ( 16. 2D /las f/ie explicit expression 



c(X) = - \\( P S C , m U u m U 2 , 



,T] m U m )\\, 

where the matrices S c and I/,- are defined as in ( I6.3D . 

Remark 6.1. from ( 16.41 ) we have the relative condition number 

(WQWfSc, myuu \\A 2 \\ F U 2 , ■ • • , \\A m \\ F U„ 



(6.4) 



CreliX) 
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6.2. The real case 

In this subsection we consider the real case, i.e., all the coefficient matrices A,, Q of Eg. ( 11. Il l 
are real. In such a case the corresponding solution X is also real. Completely similar arguments 
as Theorem l6.1| gives the following theorem. 



Theorem 6.2. Let A,, Q be real, c(X) be the condition number defined by ( 16.21 ). Then c(X) has 
the explicit expression 

1 



c(X) = - || (pS r , rjxUu mUi, ■ ■ ■ ,r]mU m ) 



where 



S, 



ZSin p(Jl r 



[(AI + Xy l Xi Af ® [(AI + Xy l X~-AX"A Pi dA 



1 + 

V i=l 

U t = S r [I ® (AjX Pl ) + ((Ajx p ') <g> /)ITJ, i = 1,2, • • ■ ,m. 
Remark 6.2. /« f/ie reaZ ca.se f/ze relative condition number is given by 

(\\Q\\ F S r , \m\pUu ||A 2 || F I/ 2 ,---,lklmllFf/ H 



c re ,(X) = 



\\X\\ f 



7. Numerical Examples 

To illustrate the results of the previous sections, in this section three simple examples are 
given, which were carried out using MATLAB 7.1. For the stopping criterion we take s^+^X) = 

m 

\\X - 2 A*XP>Ai - 1\\ < l.Oe- 10 

Example 7.1. We consider the matrix equation 

X-A\X^-Ai -A* 2 X L iA 2 = /, 

with 



+ 2 x 1(T 2 



-A, A 2 



\ + 3 x 10~ 2 



-A, A 



2 0.95 
1 



IIAII l|A|| 
Suppose that the coefficient matrices A\ and Ai are perturbed to A,- = A; + AA,-, i = 1,2, where 

10-> T 3 x lO^ 1 T 

M i = „^t . + C X AA 2 = — = — —(C T + C) 



nc r + cir 



l|C r + C|| 



and C is a random matrix generated by MATLAB function randn. 

We now consider the corresponding perturbation bounds for the solution X in Theorem U.l 
and Theorem \4.2\ 

The conditions in Theorem\4.1\are 



conl = 2(/3-s)-b>0, con2 = 8+\\AQ\\-^(l-p i )/3 p >\\A i \\ 2 >0, 
conl, = Z? 2 - 406-s)O + ||Aei|) >0. 



i=1 
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The conditions in Theorem \4.2\ are 
conA — 1 — o~ > 0, con5 



d-o-) 2 



e>0. 



By computation, we list them in Table\l\ 

Table 1: Conditions for Example l7.1l with different values of j 



j 


A 


5 


6 


7 


conl 


1.1139 


1.1141 


1.1141 


1.1141 


con2 


0.9358 


0.9358 


0.9357 


0.9357 


con3 


0.8751 


0.8755 


0.8756 


0.8756 


conA 


1.0000 


1.0000 


1.0000 


1.0000 


con5 


0.7955 


0.7957 


0.7957 


0.7957 



The results listed in Table\l\show that the conditions of Theorem \4.1\ and Theorem \4.2\ are 
satisfied. 



By Theorem \4.1\ and Theorem \4.2\ we can compute the relative perturbation bounds ft, ft = 
jpj^jj, respectively. These results averaged as the geometric mean of 10 randomly perturbed runs. 
Some results are listed in Table\2\ 

Table 2: Results for Example l7. 1 I with different values of j 



j 


4 


5 


6 


7 


ira 
ft 
ft 


3.9885 x 10~ 5 
1.9765 x 10- 4 
6.5069 x 10~ 5 


5.1141 x 10~ 6 
2.3869 x 10~ 5 
7.6524 x 10~ 6 


3.6513 x 10~ 7 
1.8133 x 10- 6 
6.0514 x 10" 7 


4.6136 x 10 s 
2.1028 x 10" 7 
6.9911 x 10~ 8 



The results listed in Table [2] show that the perturbation bound £2 given by Theorem \4.2\ is 
fairly sharp, while the bound %\ given by Theorem \4.1\ which does not depended on the exact 
solution is conservative. 



Example 7.2. We consider the matrix equation 



x - a:x {15 Ai - a;x°- 25 a 2 = i, 



with 



+ 2 x 10~ 2 



IIAII 



-A, A 2 



+ 3 x 10~ 2 
\\A\\ 





{ 2 


1 








^ 


-2 


1 


2 


1 








-A, A = 





1 


2 


1 













1 


2 


1 




, 








1 


2 , 



Choose X\ — A, X 2 — 2A. Let the approximate solution Xk of X be given with the iterative 
method ( 12.71 ), where k is the iterative number. 

The residual R{Xk) = I + A*X® ,5 Ai + A* 2 X^ 25 A2 — Xk satisfies the conditions in Theorem \5.1\ 

15 



By Theorem \5.1\ we can compute the backward error bound for X k 

\\x k -x \\< nWi%)\l 

where 

oil y mi V-l j| 1111 

P- = 1 k J ==» i = l+ll^ 1 |IIM)l|-(0.5|ft I AA _I || 2 +0.75||Z; i A 2 ^" 5 || 2 ). 

+ y0 2 -4||x- 1 ini/?(^)ii 

Some results are listed in Tabl$3\ 

Table 3: Results for Example l7.2l with different values of k 



k 


8 


10 


12 


14 


\\x k -x\\ 


6.1091 x 1(T 4 


4.0865 x 10~ 5 


2.6837 x 10~ 6 


1.7372 x 10~ 7 


iM&£)\\ 


7.2094 x 10~ 4 


4.8224 x 10~ 5 


3.1670 x 10" 6 


2.0506 x 10~ 7 



The results listed in Table\3\show that the error bound given by Theorem \5.1\ is fairly sharp. 
Example 7.3. We study the matrix equation 

X - A\X^A\ -A* 2 X i >A 2 = Q, 

with 

. ( 0.55 + 10"* \ , 1,^/11 



)' z 2 " * \ 1 

By Remark \4.2\ we can compute the relative condition number c re i(X). Some results are listed in 
Tabled 

Table 4: Results for Example l7.3l with different values of k 



k 


1 


3 


5 


7 


9 


Crel(X) 


1.1888 


1.1025 


1.1019 


1.1019 


1.1019 



The numerical results listed in the second line show that the unique positive definite solution 
X is well- conditioned. 
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